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Using the Ultra-relativistic Quantum Molecular Dynamics (UrQMD) model, we investigate az¬ 
imuthal correlations in p-Pb collisions at ^/s NN = 5.02 TeV. Comparison with the experimental 
data shows that UrQMD can not reproduce the multiplicity dependence of 2- and 4-particle cumu- 
lants, especially the transition from positive to negative values of C2{4} in high multiplicity events, 
which has been taken as experimental evidence of collectivity in p-Pb collisions. Meanwhile, UrQMD 
can not qualitatively describe the differential elliptic flow, V 2 (pt), of all charged hadrons at vari¬ 
ous multiplicity classes. These discrepancies show that the simulated hadronic p-Pb systems can 
not generate enough collective flow as observed in experiment, the associated hadron emissions are 
largely influenced by non-flow effects. However, the characteristic V 2 {pt) mass-ordering of pions, 
kaons and protons is observed in UrQMD, which is the consequence of hadronic interactions and not 
necessarily associated with strong fluid-like expansions. 

PACS numbers: 25.75.Ld, 25.75.Gz, 25.75.-q, 24.10.Lx 


I. INTRODUCTION 

The relativistic heavy ion collisions at Relativistic 
Heavy Ion Collider (RHIC) and the Large Hadron Col¬ 
lider (LHC) have provided strong evidences for the cre¬ 
ation of the Quark-Gluon Plasma (QGP) [1- ]. One 
of the crucial observables is the azimuthal anisotropy 
of the transverse momentum distribution for produced 
hadrons [5] . As a signature of the collective flow, it pro¬ 
vides important information on the Equation of State 
(EoS) and the transport properties of the QGP [6-11]. 
Usually, the anisotropy is characterised by the Fourier 
flow-coefficients [12]: 

v n = (cos[n(ip - & n )), (1) 

where ip is the azimuthal angle of the emitted hadrons, 
\E'„ is the n th -order participant (symmetry) plane angle 
and () denotes an average of all particles in all events. 
The second Fourier flow-coefficient V 2 is called elliptic 
flow, which is associated with the initial elliptic overlap 
region of the two colliding nuclei. In the past decades, 
most attention had been paid to the elliptic flow V 2 , 
which has been systematically measured and studied at 
the Super Proton Synchrotron (SPS) [ 3], RHIC [14- 
17], and the LHC [18-20] (for a summary, please also 
refer to [21, 22]). More recently, it was realised that 
the higher order flow-coefficients are equally important, 
which provide information on the fluctuating initial pro¬ 
files of the created QGP [23-37]. 
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The measurements of azimuthal correlations in 
v /s NN = 5.02 TeV p-Pb collisions at the LHC were 
originally aimed to provide reference data for the high 
energy Pb-Pb collisions, especially on the cold nu¬ 
clear matter effects. However, a large amount of un¬ 
expected collective behaviors have been discovered by 
the ALICE, ATLAS and CMS Collaborations. For in¬ 
stance, a symmetric double ridge structure on both 
near- and away-side has been observed in high multi¬ 
plicity p-Pb collisions by the ALICE Collaboration [38] . 
In addition, the CMS Collaboration has showed com¬ 
patible results between multi-particle (including 4-, 6- 
and 8-particles) and all-particle correlations with Lee- 
Yang Zero’s (LYZ) [39], which corresponds to u 2 {4} ss 
u 2 {6} ~ u 2 {8} ~ v 2 {LYZ} [40] (These results have also 
been confirmed by the ATLAS [ ] and ALICE Col¬ 
laborations [ 2]). Recently, the measurements of az¬ 
imuthal correlations have been extended to identified 
hadrons [43, 44]. A u 2 mass-ordering feature, which 
says that the differential elliptic flow at low transverse 
momentum region monotonically increases with the de¬ 
crease of hadron mass, has been observed among pions, 
kaons and protons in high multiplicity events [43]. Sim¬ 
ilarly, the CMS Collaboration found the mass-ordering 
between Kg and A(A), which showed the u 2 of Kg is 
larger than the one of A(A) at lower pp, followed by 
a crossing at pt ~ 2 GeV [44]. Many of these experi¬ 
mental measurements have been semi-quantitatively de¬ 
scribed by (3+l)-d hydrodynamic simulations from sev¬ 
eral groups [45-49], which supports the experimental 
claim that large collective flow has been developed in 
small p-Pb systems. 

In Au- Au or Pb-Pb collisions at RHIC and the LHC, 
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the collective flow mainly develops in the QGP phase 
since the QGP fireball has long enough lifetime to de¬ 
velop the momentum anisotropy until the saturation 
is almost reached [6, 50, 51]. Meanwhile, a certain 
amount of collective flow is further accumulated in the 
hadronic stage through the microscopic rescatterings, 
which leads to the v 2 mass-ordering among various 
hadron species [52-55]. Compared with Au-Au or Pb- 
Pb collisions, the smaller systems created in p-Pb col¬ 
lisions have much shorter lifetime. As a result, the mo¬ 
mentum anisotropy is not likely to reach saturation even 
if the QGP has been created. The measured azimuthal 
correlations in p-Pb collisions might be largely influ¬ 
enced by the hadronic evolution. On the other hand, 
non-flow effects (e.g. from hadron resonance decays) are 
significantly enhanced for a smaller system with much 
lower particle yields, which also contribute to 2-particle 
correlations [21]. 

With an assumption of early tliermalization for the 
created p-Pb systems, hydrodynamics simulates the 
evolution of both QGP and hadronic phases, and asso¬ 
ciate the azimuthal correlations of all charge and iden¬ 
tified hadrons with the collective expansion of the sys¬ 
tems [15-49]. In this paper, we assume that the high 
energy p-Pb collisions do not reach the threshhold of 
the QGP formation, only pure hadronic systems are 
produced. We utilize a hadron cascade model Ultra- 
relativistic Quantum Molecular Dynamics (UrQMD) [56— 
58] to simulate the evolution of the hadronic matter and 
then study the azimuthal correlations of final produced 
hadrons. Our research focuses on two aspects: (1) inves¬ 
tigating whether pure hadronic interactions could gen¬ 
erate the observed flow signatures in high multiplicity 
events; (2) studying the mass-ordering of 2-particle cor¬ 
relations in pure hadronic p-Pb systems. 

The paper is organized as follows. Section II briefly 
introduces the UrQMD model. Section III outlines the 2- 
and 4-particle Q-Cumulant method. Section IV com¬ 
pares experimental measurements with the UrQMD cal¬ 
culations on 2- and 4- particle azimuthal correlations, 
including centrality dependence and transverse momen¬ 
tum dependence. Section V summarizes and concludes 
this work. 


II. UrQMD HADRON CASCADE MODEL 

UrQMD is a microscopic transport model to describe 
hadron-hadron, hadron-nucleus and nucleus-nucleus 
collisions at relativistic energies, based on the Boltz¬ 
mann equations for various hadron species [56-58]. It 
has successfully described the soft physics at the AGS 
and SPS energies, where the created systems are domi¬ 
nated by strongly interacting hadrons. 

In UrQMD, the initial hadron productions are modeled 
via the excitation and fragmentation of strings. For 


Event class 

Nch ( \v\ < 1.0) 

0-5% 

>72 

5-10% 

60-72 

10-20% 

47-60 

20-40% 

38-47 

40-60% 

17-23 

60-100% 

<17 


TABLE I. Event class determination in UrQMD according to 
the number of all charged hadrons within |r/| < 1.0. 


Event class 

N ch (2.8 < V < 5.1) 

0-20% 

>88 

20-40% 

54-88 

40-60% 

30-54 

60-80% 

13-30 

80-100% 

< 13 


TABLE II. Event class determination in UrQMD according to 
the number of all charged hadrons within 2.8 < rj < 5.1. 


higher collision energies above = 10 GeV, the 

PYTHIA mode [59] is implemented to describe the hard 
processes and the related hadron productions. The clas¬ 
sical trajectories of the produced hadrons are then sim¬ 
ulated through solving a large set of Boltzmann equa¬ 
tions with flavour dependent cross sections. In the later 
version, UrQMD contains 55 baryon and 32 meson species 
with masses up to 2.25 GeV, supplemented by the corre¬ 
sponding antiparticles and isospin-projected states [58]. 
The elementary cross sections in the collision terms 
are either fitted from the experimental data or calcu¬ 
lated via models e.g. a modified additive quark model 
(AQM). For two closely propagating hadrons, whether 
or not a collision happens is determined by a critical dis¬ 
tance associated with the related cross section. When 
all elastic and inelastic collisions cease and all unstable 
hadrons have decayed into stable hadrons, the system 
is considered to reach kinetic freeze-out. UrQMD then 
outputs the momentum and position information of the 
final produced hadrons. 

In this paper, we implement UrQMD version 3.4 to sim¬ 
ulate the evolution of the assumed hadronic systems cre¬ 
ated in high energy p-Pb collisions. The simulations are 
executed in the equal speed system of two colliding nu¬ 
cleons with = 5.02 TeV. Correspondingly, the out¬ 
put information for final produced hadrons are defined 
in the centre-of-mass frame. In order to compare with 
the experimental data in the laboratory frame, we make 
a transformation between the centre-of-mass frame and 
the laboratory frame, which shifts the rapidity by 0.465. 
Following the related experimental papers [42, 43, 60], 
the UrQMD outputs are divided into several multiplic¬ 
ity classes, determined by the number of all charged 
hadrons N c h within a pseudorapidity range r/ < 1 or 
2.8 < < 5.1. The N c h values in these two centrality 

definitions are shown Table I and II. The pseudorapidity 
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FIG. 1. Pseudorapidity density of all charged hadrons in 
minimum bias p-Pb collisions at ^s NN = 5.02 TeV, mea¬ 
sured by ALICE [61] and calculated from UrQMD. 


density of all charged hadrons as a function of pseudo¬ 
rapidity in minimum bias p-Pb collisions is presented 
in Fig. 1. 


III. ANALYSIS METHOD AND DEFINITIONS 

In this paper, the azimuthal correlations are calcu¬ 
lated using 2- and 4-particle Q-Cumulant method [62, 
63], which were used in experiment at RHIC [64] and 
the LHC [18, 33, 40, 65]. In this method, both 2- 
and multi-particle azimuthal correlations are analyti¬ 
cally expressed in terms of a Q-vector, which is defined 
as: 

M 

Qn=Y. eiUVi ’ ( 2 ) 

i— 1 

where M is the multiplicity of the Reference Flow Parti¬ 
cles (RFPs) and p is their azimuthal angle. The single¬ 
event average 2- and 4-particle azimuthal correlations 
can be calculated via: 

,nv \Qn\ 2 -M 

{ ' M(M - 1) ’ 

\Qn\ 4 + |Q 2 »| 2 - 2 ■ Rc[q 2 „q;q*] 

W M(M — 1)(M — 2)(M — 3) [> 

0 2(M-2)- |Q„| 2 -M(M-3) 

M(M — 1)(M — 2) 


here ( ) stands for the average over all particles in a 
single event. 

The 2- and 4-particle cumulants could be achieved as: 

c„{2} = ((2)), 

c„{4} = «4))-2.<(2» 2 , 


here (( )) denotes the average over all particles over all 
events. 

In order to proceed with the calculation of the differ¬ 
ential flow of the Particles Of Interests (POIs), the p n 
and q n vectors for specific kinematic range and/or for 
specific hadron species are needed: 

m p 

i=1 ( 5 ) 

m q v ' 

i=i 

where m p is the total number of particles labeled as 
POIs, uiq is the total number of particles tagged both 
as RFP and POI. 

The single-event average differential 2- and 4-particle 
azimuthal cumulants are calculated as: 

/ 2 '\ _ PnQn ~ m q 
m.pM — m q ’ 

(4) = \pnQnQ* n Q* n ~ <hnQ* n Ql ~ PnQnQln 

- 2 • M Pn Q* n - 2 • m q \Q n \ 2 + 7 ■ q n Q* n - Q n q* n 

+ q2nQ*2n + 2 ’ PnQn + 2 ’ ™qM - 6 • TO,]. 

/[( m p M - 3 m q )(M - 1 )(M - 2 )(M - 3)] 

(6 ) 

For detectors with uniform azimuthal acceptance the 
differential 2- and 4-particle cumulants are given by: 

dn{ 2} = ((2')), 

dn{4} = ((4'))-2«2')>«2)). 


Finally the estimated differential flow from 2- 

and 4-particle correlations are given by: 


^n{2}(p T ) 


d n { 2} 

V / c42}’ 


v«{4}(pt) = 


dn{ 4} 

(-Cn{4}) 3 / 4 ' 


( 8 ) 


Unfortunately, the v n obtained from the 2-particle 
Q-Cumulant contains contributions from so-called non¬ 
flow effects, which are additional azimuthal correla¬ 
tions between the particles due to e.g. resonance de¬ 
cays, jet fragmentation, and Bose-Einstein correlations. 
They can be suppressed by appropriate kinematic cuts. 
For instance, one can introduce a pseudorapidity gap 
between the particles in the 2-particle Q-Cumulant 
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FIG. 2. m t spectra of pions, kaons and protons in p-Pb collisions at ^/s NN = 5.02 TeV measured by ALICE [60] and 
calculated from UrQMD. Here the multiplicity class determination in UrQMD is based on Table II . 


method [65]. Accordingly, the whole event is divided 
into two sub-events, A and B , which are separated by 
a |A»j| gap. This modifies Eq. (3) to: 


(2) At/ 


Qn ■ Qt 

M a ■ M b ’ 


(9) 


where Q„ and Q// are the flow vectors from sub-event 
A and B, M 4 and Mb are the corresponding multiplic¬ 
ities. 

The 2-particle Q-Cumulant with a |A?y| gap is given 
by: 


c n{2, |A 77 I} = {{2})a?7- (10) 


For the calculations of differential flow with a pseudo¬ 
rapidity gap, there is no overlap of POIs and RPs if we 
select RPs from one subevent and POIs from the other. 
This modifies Eqs. ( 6 ) to: 


(2 ) a r) 


Pn,AQn,B 
m PtA M B 1 


( 11 ) 


and we get the differential 2 -particle cumulant as: 

d n {2, lAr?]} = «2'» a „. ( 12 ) 


Finally, the differential flow from 2-particle cumulant 
can be obtained by inserting the 2 -particle reference 
flow (with 77 gap) to the differential 2 -particle cumu¬ 
lant: 


^n{2, |Ar/|}(p T ) 


<4{2,1A 77 I } 
\Jc n { 2: |A?7|} 


(13) 


In this paper, the second and third Fourier flow- 
coefficients are evaluated using above equations, by set¬ 
ting the n = 2 and 3, respectively. 


IV. RESULTS AND DISCUSSIONS 

This section mainly investigates the second and third 
Fourier flow-coefficients with cumulants in p-Pb col¬ 
lisions at = 5.02 TeV. Before studying the 2- 

particle and 4-particle correlations, it is important to 
check the single hadron information. Fig. 1 plots the 
pseudorapidity density of all charged hadrons dN c h/dr] 
in minimum bias p-Pb collisions. In general, UrQMD 
roughly describes the forward-backward asymmetry of 
the dN^/dr] curve within |? 7 | < 2. At midrapidity, 
dN c h/dr] from UrQMD are close to the ALICE measure¬ 
ments, but about 5% lower than the experimental val¬ 
ues. 

Figure 2 plots the wt spectra of pions, kaons and 
protons in high energy p-Pb collisions. It is generally 
believed that, in the absence of radial flow, tot spectra 
as a function of tot — m 0 {mo stands for the mass of the 
hadron and tot = \/Vr + m o ) satisfies the tot scaling , 
where the slope of the spectra is independent of hadron 
species [ 66 ]. Such tot scaling has been observed in p- 
p collisions at = 200 GeV [ 66 , 67]. In heavy ion 

collisions at the SPS energies and above, the tot scaling 
is broken, which provides evidence for the development 
of strong radial flow in the hot QCD systems [66-70]. 

In high energy p-Pb collisions at ,ys NN = 5.02 TeV, 
the ALICE measurements in Fig. 2 show that mx scal¬ 
ing is broken at 0-20% and 20-40% multiplicity classes 
where the measured protons spectra are flatter than 
kaons ones 1 . This provides evidence for the develop¬ 
ment of radial flow in high multiplicity events. The 


1 The pion spectra is largely influenced by resonance decays at 
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FIG. 3. (Color online) 02 ( 2 } of all charged hadrons in p-Pb collisions at ys NN = 5.02 TeV, calculated from UrQMD (Left) 
and measured by ALICE (Right) [ 2]. The circle, square and diamond markers represent various pseudorapidity gap cuts 
without r/ gap, with gap | A?/ > 0.4 and | Ar/\ > 1.0, respectively. Here the multiplicity class determination in UrQMD is based 
on Table I. 
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FIG. 4. (Color online) C2{4} of charged particles in p-Pb 
collisions at ^/s NN = 5.02 TeV, calculated from UrQMD and 
measured by ALICE [ 12] . Here the multiplicity class deter¬ 
mination in UrQMD is based on Table I. 


UrQMD calculations in Fig. 2 also present a weak bro¬ 
ken of the ttit scaling, but show steeper spectra for pi- 
ons, kaons and protons when compared with the ALICE 
curves. This indicates that the assumed hadronic p-Pb 
systems could not accumulate sufficient radial flow as 
observed in experiment. 


lower raj — mo, which break the pion’s rax scaling even for the 
case without radial flow 


With brief investigations of the single hadron data, 
we now focus on studying azimuthal correlations in 
high energy p-Pb collisions. Figure 3 presents the cen¬ 
trality dependence of the 2 -particle cumulant of the 
second Fourier flow-coefficient C 2 { 2 }, calculated from 
UrQMD hadron cascade model (left) and measured by 
the ALICE collaboration (right). For various pseudora¬ 
pidity gaps, c 2 {2} from UrQMD exhibit decreasing trend 
from peripheral (low multiplicity events) to central col¬ 
lisions (high multiplicity events), which agrees with the 
expectation of the azimuthal correlations not associated 
with the symmetry plane, the so-called non-flow effects. 
As the pseudorapidity gap increases, the magnitudes of 
c 2 {2} become weaker for both ALICE and UrQMD, which 
illustrates that non-flow effects, usually few-particle cor¬ 
relations from resonance decays and jets, are suppressed 
by a large pseudorapidity gap. When the pseudorapid¬ 
ity gap | A 77 I is larger than 1.0, c 2 {2} from ALICE show 
much weaker centrality dependence, which is suggested 
as one of the hints for collective expansion in the cre¬ 
ated p-Pb systems. However, c 2 {2} from UrQMD still 
present strong centrality dependence for |A?y| > 1.0, 
showing a typical non-flow behavior. Usually, the non¬ 
flow effects between 2 -particle correlations, denoted as 
S n , behave as S n ~ 1 /M where M is the multiplicity. 
The decreasing trend of c 2 {2} with the increase of mul¬ 
tiplicity indicates that UrQMD hadronic expansion could 
not generate enough flow in a small p-Pb system, non¬ 
flow effects are still pretty large even for the case with 
a large pseudorapidity gap cut |Ar/| > 1.0. 

To better understand the hadronic systems simulated 
by UrQMD, we investigate the 4-particle cumulant of the 
second Fourier flow-coefficient c 2 {4}, which is equal to 
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FIG. 5. (Color online) 03 ( 2 } of all charged hadrons in p-Pb collisions at ^/s NN = 5.02 TeV, calculated from UrQMD(Left) and 
measured by ALICE (Right) [42]. Here the multiplicity class determination in UrQMD is based on Table I. 
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—r> 2 { 4} 4 and expected to be less sensitive to non-flow ef¬ 
fects. Figure 4 plots the centrality dependence of C 2 { 4 } 
of all charged hadrons in p-Pb collisions at y/s^ = 
5.02 TeV. Both the UrQMD and ALICE results show 
that c 2 {4} increase with the decrease of multiplicity 
from semi-central to peripheral collisions. For the most 
central collisions (<10%), C 2 { 4 } from ALICE exhibits 
a transition from positive to negative values, indicat¬ 
ing the creation of flow-dominated systems in the high 
multiplicity events. However, C 2 { 4 } from UrQMD keeps 
positive for all available multiplicity classes, including 
the most central collisions. As a result, real values of 
u 2 {4} can not be extracted in UrQMD for all centrality 
bins. This comparison further illustrates the difference 
between the p-Pb systems created in experiment and 
simulated by UrQMD. The hadron emissions from UrQMD 
are largely influenced by non-flow effects. Without the 
contributions from the initial stage and/or the QGP 
phase, the measured flow-like 4-particle correlations in 
high multiplicity events can not be reproduced by a mi¬ 
croscopic transport model with only hadronic scatter¬ 
ings and decays. 

In Figure 5, we further study the 2-particle azimuthal 
correlations for the third Fourier flow-coefficient C 3 { 2 }. 
The UrQMD calculations and the ALICE measurements 
with various pseudorapidity gaps are respectively shown 
in the left and right panels of Fig. 5. Similar to c 2 {2} 
in Fig. 3, C 3 { 2 } also decreases with the increase of |Ar/|. 
For the ALICE measurement, 03 ( 2 } keeps positive for 
all pseudorapidity gaps, which leads to real values of 
triangular flow v 3 {2} (d 3 {2} = ^/c 3 { 2 }) as measured 
in [42]. Considered that non-flow effects are largely 
suppressed by a large pseudorapidity gap, the measured 
c 3 {2} at |Ary| > 1.0 (and the associated triangular flow 


u 3 {2}) is possibly mainly caused by collective expan¬ 
sion and reflects initial state fluctuations of the p-Pb 
systems. In contrast, c 3 {2} from UrQMD turns to neg¬ 
ative for |Ar/| > 0.4 and |A 77 I > 1.0, which does not 
produce a real value of V 3 { 2 } 2 . The fact that UrQMD 
could not generate the experimentally observed trian¬ 
gular flow, together with the results shown Figs. 3-5, 
strongly indicates that p-Pb systems from UrQMD con¬ 
tain large non-flow azimuthal correlations. 

Following Eq. (13), we calculate the second Fourier 
flow-coefficient as a function of transverse momentum, 
V2 (pt ) 1 for the UrQMD simulations at multiplicity class 
0-20%, 20-40%, 40-60% and 60-100%. Figure 6 shows 
that u 2 (pt) monotonically increases from high to low 
multiplicity class, which agrees with the trend of c 2 { 2 } 
shown in Fig. 3 (c 2 {2} is the square of integrated V2{2}). 
Meanwhile, u 2 (pt) from UrQMD increases with the in¬ 
crease of p T and show strong sensitivity to the pseu¬ 
dorapidity gap. The observed large pseudorapidity gap 
suppression of u 2 (pt) indicates that non-flow effects are 
large in UrQMD, as already shown in Fig. 3-5. 

Figure 6 also shows that UrQMD can not correctly re¬ 
produce the shape of the experimental V 2 (pt) curves, 
when implemented the same pseudorapidity gap | A 77 I > 
1.0. It underpredicts the data at lower px and overes¬ 
timates the data above 1 GeV. Compared with the in¬ 
tegrated i> 2 , the differential elliptic flow u 2 (pt) contains 
more information on the evolving system, which reflects 
the interplay between radial and elliptic flow. The tot 
spectra in Fig. 2 has already shown that UrQMD can not 


2 We also find that 03 ( 4 } only shows positive values, just as 
C2{4}, which does not produce a real value of 1 ) 3 ( 4 }. 
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FIG. 6. (Color online) V 2 (pt) of all change hadrons in p-Pb collisions at ^s NN = 5.02 TeV, calculated from UrQMD and 
measured by ALICE [ 13]. Here the multiplicity class determination in UrQMD is based on Table II. 


produce sufficient radial flow as observed in experiment. 
The insufficient radial flow, together with the insuffi¬ 
cient flow anisotropy accumulation (shown in Fig. 3-5) 
leads to the fact that UrQMD could not reproduced the 
v- 2 (pr) curves measured by ALICE. 

Figure 7 investigates azimuthal correlations of iden¬ 
tified hadrons in high energy p-Pb collisions. The right 
panels present the ALICE measurements with two dif¬ 
ferent multiplicity classes [43, 44], which show a char¬ 
acteristic feature of V 2 (pt) mass-ordering among pi- 
ons, kaons and protons. In the past research, hydro- 
dynamic simulations from several groups have system¬ 
atically studied the flow data, which reproduced the 
v -2 mass-ordering feature of the p-Pb systems [46, 49]. 
In the hydrodynamic language, the radial flow further 
accumulated in hadronic stage tends to push heavier 
hadrons from lower px to higher px, leading to an en¬ 
hanced V 2 splitting between pions and protons [53, 55]. 
The observation of mass-ordering is thus generally 
believed as a strong evidence for the collective expan¬ 
sion of the p-Pb systems created in y/s^ =5.02 TeV 
collisions. 

However, the left panels of Fig. 7 shows that UrQMD 
also generate a mass-ordering for the 2-particle corre¬ 


lations among pions, kaons and protons 3 . Such mass¬ 
ordering pattern, caused by pure hadronic interactions, 
qualitatively agrees with the ones from the ALICE 
measurement [43] and from the hydrodynamic calcula¬ 
tions [ 16, 49]. In UrQMD, the unknown cross sections are 
calculated by the additive quark model (AQM) through 
counting the number of constituent quarks within two 
colliding hadrons. As a result, the main meson-baryon 
(M-B) cross sections from AQM are about 50% larger 
than the meson-meson (M-M) cross sections, leading to 
the i >2 splitting between mesons and baryons after the 
evolution of hadronic matter. Comparison simulations 
in appendix A (Fig. 9) also show that, with the M-B 
and M-M interaction channels closed in UrQMD, the V 2 
mass-ordering almost disappears. The combined results 
in Fig. 7 and 9 illustrate that the hadronic interactions 
could lead to a mass-ordering in 2-particle correlations 
among pions, kaons and protons, even for small p-Pb 
systems without enough flow generation. 


3 Due to limited statistics, we apply |A»j| > 0.2 in our calcula¬ 
tions, rather than |Ar;| > 0.8 as used in experiment. In fact, 
V2{2; A// > 0.8} from our current UrQMD simulations has large 
error bars, especially for protons. However, a tendency of V 2 
mass ordering among pions, kaons and protons is still observed. 
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FIG. 7. (Color online) V 2 (pt) of pions, kaons and protons in p-Pb collisions at ^/s NN = 5.02 TeV, calculated from UrQMD(left 
panels) and measured by ALICE(right panels) [43]. Here the multiplicity class determination in UrQMD is based on Table II. 


V. SUMMARY 

Using UrQMD hadron cascade model, we studied az¬ 
imuthal correlations in p-Pb collisions at ^/s NN = 5.02 
TeV. Comparisons with the experimental data showed 
that the p-Pb systems created in experiment are not the 
trivial hadronic systems described by UrQMD. Here, we 
summarize the main results as the following: 

(1) With large pseudorapidity gaps (|Ar/| > 1.0), 
the measured 2-particle cumulant of the second Fourier 
flow-coefficient C2{2} from ALICE shows a weak cen¬ 
trality dependence from central to semi-peripheral col¬ 
lisions. In contrast, the UrQMD calculations still present 
a strong centrality dependence for c 2 {2, |Ar/| > 1.0}. 

(2) In the most central collisions, c 2 {4} from ALICE 
exhibits a transition from positive to negative values, 
which indicates the development of strong collective 
flow in high multiplicity events. However, c 2 {4} from 
UrQMD keeps positive for all multiplicity classes, which 
does not produce u 2 {4} with a real value. 

(3) For large pseudorapidity gaps, C3{2} from UrQMD 


turns to negative values, which can not produce the 
triangular flow as observed in experiments. 

(4) UrQMD can not fit the differential flow u 2 (pt) from 
ALICE at various multiplicity classes. 

More specifically, the related experimental data of az¬ 
imuthal correlations have accumulated strong evidence 
for the development of strong collective flow in high 
multiplicity events. With the assumption that high en¬ 
ergy p-Pb collisions do not reach the threshold for the 
QGP formation and only produce trivial hadronic sys¬ 
tems, we did hadron transport simulations with UrQMD. 
We found that hadronic interactions alone could not 
generate sufficient collective flow as observed in ex¬ 
periment. Non-flow effects, e.g. from resonance de¬ 
cays and/or jet-like fragmentations, largely influence 
the hadron emissions of the UrQMD systems. In order 
to fit the measured azimuthal correlations of all charged 
hadrons in p-Pb collisions at = 5.02 TeV, the con¬ 
tributions from the initial stage and/or the QGP phase 
can not be neglected. 

In addition, we extended our study of azimuthal cor- 
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FIG. 8. rriT spectra of of pions, kaons and protons in p-Pb collisions at ys NN = 5.02 TeV, calculated from UrQMD with and 
without M-M and M-B collisions. Here the multiplicity class determination in UrQMD is based on Table II. 


relations to identified hadrons. The calculations of 
the 2-particle correlations for pions, kaons and protons 
showed that UrQMD can generate a v -2 mass-ordering 
with the characteristic feature similar to the ALICE 
measurements. Comparison runs from UrQMD with 
main hadronic scatterings turned on and off showed 
that the v -2 mass-ordering in UrQMD is mainly caused 
by hadronic interactions. The V 2 mass-ordering alone is 
not necessarily a flow signature associated with strong 
fluid-like expansions. 
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Appendix A: UrQMD comparison runs 
with/without M — M and M B collisions 

This appendix explores how hadronic interactions 
in UrQMD influence spectra and azimuthal correlations 
of identified hadrons for the hadronic p-Pb systems. 
In UrQMD, hadronic scatterings include Meson-Meson 


(M-M) collisions, Meson-Baryon (M-B) collisions and 
Baryon-Baryon (B-B) collisions. When switching off all 
of these collision channels, UrQMD simulations, in prin¬ 
ciple, consist of initial hadron productions and the suc¬ 
ceeding resonance decays, which are mainly influenced 
by non-flow effects. However, not all these collision 
channels in the current version of UrQMD (v3.4) can be 
simultaneously turned off. With B-B collision channels 
turned off, all the secondary proton-nucleon collisions 
from the initial p-Pb collisions are automatically turned 
off without proceeding any further hadron productions 
and decays. Considering the probability of B-B colli¬ 
sions is much lower than M-M and M-B collisions, we 
only turn off the M-M and M-B interaction channels for 
the UrQMD comparison runs in this appendix. 

Figure 8 plots the mx spectra of pions, kaons and 
protons in p-Pb collisions at y/s NN = 5.02 TeV, cal¬ 
culated from UrQMD with and without M-M and M-B 
interactions. In Sec. IV, we have already showed that, 
although the m-x scaling is weakly broken in UrQMD, pure 
hadronic interactions can not generate sufficient radial 
flow as observed in experiment. Here, Fig. 8 shows that 
M-M and M-B collisions only slightly change the slope 
of the m t spectra 4 . The slight broken of the mx scal¬ 
ing in UrQMD are very possibly caused by mechanisms of 
the initial hadron productions 5 . 


4 For 2.76 A TeV Pb+Pb collisions at 70-80% centrality (where 
d.N r h /dtj at mid-rapidity is about 50 [ 5] and close to dN^/drj 
in high multiplicity p-Pb collisions), hybrid model simulations 
also show that hadronic scatterings almost do not change the 
mx spectra of identified hadrons [71]. 

5 In the UrQMD simulations for p-Pb collisions at ^/s NN = 5.02 
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Compared with the ?nx spectra, the V 2 mass-ordering 
are more sensitive to hadronic interactions. For typi¬ 
cal flow-dominated systems created in high-energy An 
Au or Pb-Pb collisions, hybrid model simulations have 
shown that hadronic rescatterings dramatically increase 
the V 2 splitting between pions and protons, but only 
slightly change the itit spectra [54]. Figure 9 presents 
i> 2 {pt) of identified hadrons in high energy p-Pb col¬ 
lisions, based on UrQMD simulations in the senarios of 
with (left panels) and without (right panels) M-M and 
M-B collisions. In the cases that M-M and M-B colli¬ 


sions are turned off, the V 2 mass-ordering among pions, 
kaons and protons almost disappears, when compared 
with the cases with M-M and M-B interactions. 

In Sec. IV, detailed study of the 2-particle and 4- 
particle correlations has already shown that UrQMD 
could not generate sufficient flow as observed in experi¬ 
ments, its final hadron emissions are largely influenced 
by non-flow effects. The comparisons runs in Fig. 9 il¬ 
lustrate that the V 2 mass-ordering can be explained as 
the consequence of hadronic interactions, which is not 
necessarily associated with strong fluid-like expansions. 
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